Demonstration of Two-Qubit Algorithms with a Superconducting Quantum Processor 
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By harnessing the superposition and entan- 
glement of physical states, quantum comput- 
ers could outperform their classical counter- 
parts in solving problems of technological impact, 
such as factoring large numbers and searching 
databases^'^. A quantum processor executes al- 
gorithms by applying a programmable sequence 
of gates to an initialized register of qubits, which 
coherently evolves into a final state containing 
the result of the computation. Simultaneously 
meeting the conflicting requirements of long co- 
herence, state preparation, universal gate oper- 
ations, and qubit readout makes building quan- 
tum processors challenging. Few-qubit proces- 
sors have already been shown in nuclear magnetic 
resonanceSii^i^i^, cold ion trap'"'* and optical^ sys- 
tems, but a solid-state realization has remained 
an outstanding challenge. Here we demon- 
strate a two-qubit superconducting processor and 
the implementation of the Grover search^" and 
Deutsch Jozsa-*^-*^ quantum algorithms. We em- 
ploy a novel two-qubit interaction, tunable in 
strength by two orders of magnitude on nanosec- 
ond time scales, which is mediated by a cavity 
bus in a circuit quantum electrodynamics (cQED) 
architecture-*^^'-*^^. This interaction allows genera- 
tion of highly-entangled states with concurrence 
up to 94%. Although this processor constitutes 
an important step in quantum computing with 
integrated circuits, continuing efforts to increase 
qubit coherence times, gate performance and reg- 
ister size will be required to fulfill the promise of 
a scalable technology. 

Over the last decade, superconducting circuitsii have 
made considerable progress on all the requirements nec- 
essary for an electrically-controlled, solid-state quantum 
computer. Coherence time a^^i^^ have risen by three 
orders of magnitude to '--^ 1 fj,s, single-qubit gatesi^iii 
have reached error rates of 1%, engineered interac- 
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have produced two-qubit entanglement at 
a level of 60% concurrence^^, and qubit rcadou t^^'^'^i^^ 
has attained measurement fidelities ~ 90%. However, 
combining these achievements in a single device remains 
challenging. One approach to integration is the quantum 
bus architectur o^^'^^1^^ , which uses an on-chip transmis- 
sion line cavity to couple, control, and measure qubits. 



We augment the architecture in Ref. [261 with flux-bias 
lines that tune individual qubit frequencies, permitting 
single-qubit phase gates. By pulsing the qubit frequen- 
cies to an avoided crossing where a CTz (8" ctz interaction 
turns on, we are able to realize a two-qubit conditional 
phase (c-Phase) gate. Operation in the strong-dispersive 
regime^^ of cQED allows joint readout^^ that can effi- 
ciently detect two-qubit correlations. Combined with 
single-qubit rotations, this enables tomography of the 
two-qubit state. Through an improved understanding of 
spontaneous emission^S and careful microwave engineer- 
ing, we are now able to combine state-of-the-art 1 fis 
coherence times into a two-qubit device. This allows suf- 
ficient time to concatenate ~ 10 gates, realizing simple 
algorithms with fidelity greater than 80%. 

Our processor, shown in Fig. [T^, is a 4-port super- 
conducting device comprising two transmon qubit a^^i'^° 
(Ql and Qr) inside a microwave cavity bus, and flux- 
bias lines proximal to each qubit. The cavity, nor- 
mally off-resonance with the qubit transition frequen- 
cies /l and /r, couples the qubits by virtual photon 
exchange and shields them from the electromagnetic 
continuum. As previously demonstrated^^, microwave 
pulses resonant with /l or /r applied to the cavity input 
port provide frequency-multiplexed single-qubit x- and 
?/-rotations with high fidclity^^ and selectivity. Pulsed 
measurement of the homodyne voltage Vh on the output 
port of the cavity provides qubit readout. The remaining 
two ports create local magnetic fields that tune the qubit 
transition frequencies. Each qubit has a split Josephson 
junction, so its frequency depends on the flux <f> through 
the loop according to hf « y/8E^Ef^''\coa{'i>/^o)\-Ec, 
where Ec is the charging energy, Ej^'^^ is the maximum 
Josephson energy, h is Planck's constant, and $o is the 
flux quantum. By employing short-circuited transmis- 
sion lines with a bandwidth from dc to 2 GHz, we can 
tune /l and /r by many GHz using room temperature 
voltages Vl and Vr . Static tuning of qubit transitions 
using the flux-bias lines is demonstrated in Fig. [T)3. 

The spectrum of single excitations (Fig. [TJs) shows the 
essential features of the cavity-coupled two-qubit Hamil- 
tonian and allows a determination of the relevant system 
parameters (see Methods). When the qubits are tuned to 
their maximum frequencies, point I, they are far detuned 
from the cavity and from each other, so that interactions 
are small. This point is therefore used for state prepa- 
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Figure 1: Two-qubit cQED device, and cavity/qubit characterization, a, Optical micrograph of 4-port device with 
a coplanar waveguide cavity bus coupling two transrnon qubits (insets), and local flux-bias lines providing fast qubit tuning. 
Microwave pulses at the qubit transition frequencies /l and /r drive single-qubit rotations, and a pulsed measurement of 
the cavity homodyne voltage Vh (at frequency /c) provides two-qubit readout. The flux-bias lines (bottom-left and top-right 
ports) are coplanar waveguides with short-circuit termination next to their target qubit. The termination geometry allows 
current on the line to couple flux through the split junctions (b, inset), b, Grey scale images of cavity transmission and of 
qubit spectroscopy as a function of Vr, showing local tuning of Qr across the avoided crossing with Ql (point III) and across 
the vacuum Rabi splitting with the cavity (point IV). Semi-transparent lines are theoretical best fits obtained from numerical 
diagonalization of a generalized Tavis-Cummings Hamiltonian''^ . Points I and II are the operating points of the processor. 
Preparation, singlc-qubit operations and measurements are performed at point I, and a c-Phase gate is achieved by pulsing 
into point II. 



ration, single-qubit rotations and measurement, in the 
computational basis |0, 0), |0, 1), |1,0), and where 
\l,r) denotes excitation level / (r) for Ql (Qr)- Oper- 
ation at this point is also desirable because it is a flux 
sweet spolji^ for both qubits, providing long coherence, 
with relaxation and dephasing times Ti l(i^) = 1.3(0.8) /zs 
and ?^2*L(R) ~ 1.8(1.2) /is, respectively. Tuning Qr into 
resonance with the cavity, point IV, reveals a vacuum 
Rabi splitting^^ from which the qubit-cavity interaction 
strength is extracted. Tuning into resonance with 
Ql, point III, shows an avoided crossing resulting from a 
cavity-mediated, qubit-qubit transverse interaction^i^''^?^ 
investigated previously^. In this work, we perform two- 
qubit gates at point II, where no interactions arc immedi- 
ately apparent on examining the onc-cxcitation manifold. 

However, a useful two-qubit interaction is revealed in 
the two-excitation spectrum, shown in Fig. As ^r, is 
swept away from point I, the non-computational higher- 
level transmon excitation |0,2) decreases more rapidly 
than the computational state |1,1), and these states 
would become degenerate at point II. But as shown 
in Fig. there is a large (160 MHz) cavity-mediated 
interaction between these levels, producing a frequency 
shift C/27r of the lower branch with respect to the sum 
/l + /r, in good agreement with a numerical diagonaliza- 
tion of the generalized Tavis-Cummings Hamiltoniaii^i 



(see Methods). 

This shift is the mechanism at the heart of our condi- 
tional phase gate. Flux pulses, adiabatic with respect to 
the |1, 1) ^ |0,2) avoided crossing, produce phase gates 
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in the computational Hilbert space. Here, (^ir = 
27r J dfir(t)dt is the dynamical phase acquired by 
and Sfir is the deviation of fir from its value at point I. 
A Vr pulse into point II such that / C('^) dt — {2n + l)7r 
with integer n implements a c-Phase, because (pn — 0oi + 
010^/ C(^) di. This method of realizing a c-Phase by adi- 
abatically using the avoided crossing between computa- 
tional and non-computational states is generally applica- 
ble to any qubit implementation with finite anharmonic- 
ity, such as transmons^^ or phase qubits^^. A similar 
approach involving higher excitation levels but with non- 
adiabatic pulses was previously proposed'^-'. The negative 
anharmonicity permits the phase gate at point II to occur 
before the onset of transverse coupling at point III. 

Control of C, by two orders of magnitude provides an 
excellent on-off ratio for the c-Phase gate. As shown 
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Figure 2: Origin and characterization of the controlled- 
phase gate, a, Flux dependence of transition frequen- 
cies from tlie ground state |0, 0) to tiie two-excitation man- 
ifold. Spectroscopy measurements (points) show an avoided 
crossing between the computational state |1, 1) and the non- 
computational state |0, 2) at point II, in good agreement 
with numerical diagonalization of the Hamiltonian (dashed 
curves), b, This avoided crossing causes the transition fre- 
quency to |1,1) to deviate from the sum of the transition 
frequencies to |0, 1) and |1,0). c, The coupling strength 
^/27r = /oi -I- /lo — /ii of the effective (g) interaction, 
obtained both from spectroscopy (solid curve) and from time- 
domain experiments (points) (see text for details). Numerical 
diagonalization and perturbation theory (Supplementary In- 
formation) for 3-level transmons agree reasonably with data. 
The perturbation calculation diverges at the avoided cross- 
ing. Perturbation theory for 2-level qubits gives the wrong 
magnitude and sign for (, and demonstrates that the higher 
transmon excitations are necessary for the interaction. Time- 
domain measurement and theory both give C/27r — 1.2 MHz 
at point I. The tunability of ( over two orders of magnitude 
provides an excellent on-off ratio for the c-Phase gate. 



in Fig. [2j;, measurements of ( obtained from spec- 
troscopy and from time-domain experiments show very 
good agreement. The time-domain method measures 
the difference in the precession frequency of Ql in two 
Ramsey-style experiments where a VR-pulse of varying 
duration (0-100 ns) is inserted between tt/2 rotations of 
Ql, with Qr either in the ground state |0) or excited into 
state |1). Using the time-domain approach, we measure 
a residual C/^tt ~ 1.2 MHz at point I (star). The the- 
oretical C obtained by numerical diagonalization shows 
reasonable agreement with the data, except for a scale 
factor that is likely due to higher modes of the cavity, 
not included in the calculation. 

The controlled phase interaction allows universal two- 
qubit gates. As an example, we produce high-fidelity 
entangled states on demand (Fig. [3]). The pulse sequence 
in Fig.lHK generates any of the four Bell states. 
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depending on the choice of c-Phase gate cUij applied 
{cU.ij\l,r) = (— l)'''"'f |Z, r)). These gates are realized 
through fine control of the dynamical phases (j)oi and (j)io 
in a 30 ns ^R-pulse close to point II and back. We tune 
(/>oi over 2tt by making small adjustments to the rising 
and falling edges of the Vpi-pulse, and 0io with a simul- 
taneous weak VL-pulse. 

To detect the entanglement, we first reconstruct the 
two-qubit density matrix p by quantum state tomogra- 
phy using joint dispersive readou t^^i^^'^^ . A pulsed mea- 
surement of the cavity homodyne voltage Vh measures 
the operator 

M = Pia^ + I32af + /S^a^ (g> af, 

where the a are two-qubit Pauli operators^. Operation in 
the strong-dispersive regime^ii^ makes IPu] ~ \f3i\, \P2\, 
enhancing sensitivity to two-qubit correlations. A com- 
plete set of 15 linearly independent measurement opera- 
tors is built using single-qubit rotations prior to measur- 
ing M. An ensemble average of each operator is obtained 
by executing the sequence in Fig. [3^ 450,000 times. The 
15 measured values are then input to a maximum likeli- 
hood estimator— of p (see Supplementary Information). 

The inferred density matrices pmi reveal highly- 
entangled states in all four cases (Fig. [5jD-e) . We quantify 
performance using the metrics of purity, P{p) = Tr(p^), 
fidelity to the target state \tp), F{p,4') = £^nd 
concurrence^, C, computable from pmi. Note that there 
are several common definitions of fidelity in the litera- 
ture, and our definition is the square of the fidelity used 
in Refs. 22: and 28;. Values for P, F and C for the four 
cases are given in the caption to Fig. [31 These values 
significantly extend the current state of the art for solid- 
state entanglement^, and provide evidence that we have 
a high-fidelity universal set of two-qubit gates. 

One- and two-qubit gates can be concatenated to real- 
ize simple algorithms, such as Grover's quantum searcbi^ 
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Figure 3: Entanglement on demand, a, Gate sequence 
generating two-qubit entanglement and detection via quan- 
tum state tomography. Starting from |0,0), simultaneous 
7r/2 rotations on both qubits create an equal superposi- 
tion of the four computational states. A c-Phase cUij then 
phase shifts in the superposition and produces entan- 

glement. A final 7r/2 rotation on Ql evolves the entan- 
gled state into one of the four Bell states depending on the 
cUij applied, b— e, Real part of maximum-likelihood den- 
sity matrix pmi of the entangler output for c(7io, cUoo, cUn, 
and cUoi, respectively (imaginary elements of pmi are less 
than 0.03, 0.02, 0.07, 0.08). Extracted metrics for the four 
entangler outputs include purity P = 0.87 ± 0.02, 0.92 ± 
0.02, 0.88 ± 0.02, 0.79 ± 0.03, fidelity to the ideal Bell state 
F = 0.91 ± 0.01, 0.94 ± 0.01, 0.90 ± 0.01, 0.87 ± 0.02 and con- 
currence C = 0.88 ± 0.02, 0.94 ± 0.01, 0.86 ± 0.02, 0.81 ± 0.04. 
The uncertainties correspond to the standard deviation in 16 
repetitions of generation-tomography for each entangler. 



shown in Fig. U) Given a function f{x) on the set 
X G {0,1,2,3} such that f{x) — 1 except at some xq, 
where /(a;o) = — 1, this well-known algorithm can deter- 



mine xo — 21+ j with a single call of an oracle O = cU, 
which encodes f{x) in a quantum phase. 

We can examine the functioning of the algorithm 
by interrupting it after each step and performing state 
tomography. Figure S^-g clearly shows all the features 
of a quantum processor, namely the use of maximally su- 



perposed states to exploit quantum parallelism (Fig. |4);), 
the encoding of information in the entanglement between 
qubits (Fig. HJi, e), and the interference producing an 
answer represented in a final computational basis state. 
The fidelity of the final state (Fig. |4^) to the expected 
output (|1,0) for the case O = cC/io shown) is 85%. 
Similar performance is obtained for the other three 
oracles (Table HI. 

We have also programmed and executed the Deutsch- 
Jozsa algorithn>iii^. The two-qubit version of this al- 
gorithm determines whether an unknown function fi (x) , 
mapping a one-bit input to a one-bit output, is constant 
{fo{x) = or fi{x) = 1) or balanced {f2{x) = x or 
h{x) — 1 — 2;), doing so with a single call of the function. 
The algorithm applies the function once to a superposi- 
tion of the two possible inputs and employs the concept 
of quantum phase kick-back^ to encode the result in the 
final state of one qubit (here, Ql) while leaving the other 
untouched (Qr). The gate sequence realizing the algo- 
rithm and the output tomographs for the four cases are 
shown in Supplementary Fig. SI. 

The performance of both algorithms is summarized 
in Table [H Although there are undoubtedly signif- 
icant systematic errors remaining, the overall fidelity 
is nonetheless similar to that expected from the ratio 
(^ 100 ns/1 /is) of the total duration of gate sequences to 
the qubit coherence times. 

In summary, we have demonstrated the experimental 
realization of two-qubit quantum algorithms using a su- 
perconducting circuit. The incorporation of local flux 
control and joint-dispersive readout into cQED, together 
with a ten-fold increase in qubit coherence over previ- 
ous two-qubit devices, has enabled on-demand generation 
and detection of entanglement and the implementation 
of the Grover and Deutsch-Jozsa algorithms. Supercon- 
ducting circuits could eventually perform more complex 
quantum algorithms on many qubits, provided that co- 
herence lifetimes and the resulting gate fldelities can be 
further improved. 



I. METHODS 

A. Device fabrication 

A 180 nm film of Nb was dc-magnetron sputtered on 
the epi-polished surface of an R-plane corundum (a- 
AI2O3) wafer (2" diameter, 430 /xm thickness). Copla- 
nar waveguide structures (cavity and fiux-bias lines) were 
patterned by optical lithography and fluorine-based reac- 
tive ion etching of Nb. Transmon features (interdigitated 
capacitors and split junctions) were patterned on individ- 
ual 2 mm X 7 mm chips using electron-beam lithography, 
double angle evaporation of Al (20/90 nm) with interme- 
diate oxidation (15% O2 in Ar at 15 Torr for 12 min), and 
lift-off. 

A completed device was cooled to 13 mK in a '^He-'*He 
dilution refrigerator. A diagram of the refrigerator wiring 
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Figure 4: Implementation of Grover's search algo- 
rithm, a, Concatenation of single-qubit and c-Phase gates 
implementing one iteration of Grover searching. Without loss 
of generality, we have replaced the Walsh-Hadamard trans- 
formations W = RxRy^^ in the usual description of the al- 
gorithntii^ with Ry^^ rotations in order to eliminate 6 single- 
qubit rotations and complete the sequence in 104 ns. (Supple- 
mentary Fig. S3 shows the microwave and flux pulses imple- 
menting the sequence.) The orange box is the oracle O — cUij 
that encodes the solution xq = 2i + j to the search problem 
in a quantum phase. Note that the first half of the algorithm 
is identical to the entangling sequence in Fig. O while the 
second half is essentially its mirror image, b— g, Real part 
of pmi obtained by state tomography after each step of the 
algorithm with oracle O = cUio- Starting from |0,0) (b), the 
qubits are simultaneously rotated into a maximal superposi- 
tion state (c). The oracle then marks the solution, |1,0), by 
inverting its phase (d). The Ry^'^ rotation on Ql turns the 
state into the Bell state 1^'"''), demonstrating that the state 
is highly entangled at this stage. The Ry^'^ rotation on Qr, 
produces a state identical to (d) (data not shown). The appli- 
cation of cf/oo undoes the entanglement, producing a maximal 
superposition state (f). The final rotations yield an output 
state (g) with fidelity F = 85% to the correct answer, |1, 0). 



is shown in Supplementary Fig. S2. Careful microwave 
engineering of the sample holder and on-chip wirebond- 
ing across ground planes were required to suppress spuri- 



ous resonance modes on- and off-chip. Simulations using 
Sonnet® software provided guidance with this iterative 
process. The sample was enclosed in two layers of Cryop- 
crm magnetic shielding, allowing high-fidelity operation 
of the processor during unattended overnight runs. 

B. cQED Theory 

The Tavis-Cummings^'^ Hamiltonian generalized to 
multi-level transmon qubits^ is 



H = ujco-^o, + (1 

N N 

E (E^o,b%oi, + (a+«t)E 9%b)M 

ge{L,R} 3=0 j,k=0 



Here, ujc is the bare cavity frequency, 
ujQj{Ecq, E,jq) is the transition frequency for qubit q from 
ground to excited state j, and g^f. = gqnjk{Ecq, Ejq), 
with Qq a bare qubit-cavity coupling and rijk a level- 
dependent coupling matrix element. The dependence 
of these parameters on qubit charging energy Ecq and 
Josephson energy E}q is indicated. The flux control en- 
ters through = Ef^'^\cos{-K<^q/<^Q)\, with the flux 
through the qubit loop, and a linear flux- voltage relation 
<^q — aqi,VL + aqRVR + ^q.o, accouutiug for crosstalk (~ 
30%) and offsets. The above parameters are tightly con- 
strained by the combination of spectroscopy and trans- 
mission data shown (Figs. [T|d, and [2b) and transmis- 
sion data (not shown) for the QL-cavity vacuum Rabi 
splitting. By simultaneously fitting the spectra given by 
numerical diagonalization of the Hamiltonian (truncated 
to = 5 qubit levels and 5 cavity photons) to these data, 
we obtain E'-p^^^-^/h = 28.48 (42.34) GHz, ScL(R)//i = 
317 (297) MHz, 5L(R)/27r = 199 (183) MHz, and luc/2tt = 
6.902 GHz. The cavity linewidth is K/27r = 1 MHz. 
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I. SUPPLEMENTARY INFORMATION 

A. Perturbation Theory 

To gain additional insight on the large on-off ratio of 
the frequency shift Q, observed both in experiment and 
numerical diagonalization of the Hamiltonian (Fig. 2), 
we perform a perturbative analysis in the rotating-wave 
approximation. Truncating Eq. (1) at three transmon 
excitations and assuming ni2 ~ \/2 (vahd for Ej/Eq ^ 
1), we obtain the fourth-order result 
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Here, A„ = tjg^ 
,R 



^01 



uj^-^ - ujc, Si = uj^i - w^2j and 62 = uj^ 
wf'2. This expression diverges as the 0^1 transition of 
one transmon lines up with the 1^2 transition of the 
other. Assuming instead two-level qubits, the expression 
simplifies to 
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Both perturbative expressions are compared with numer- 
ical diagonalization of the Hamiltonian in Fig. 2c. The 
three-level expression shows reasonable agreement away 
from the divergence, while the two-level expression is in- 
correct in both magnitude and sign. 



B. State Tomography 

The goal of quantum state tomography is to estimate 
the density matrix p describing a quantum mechanical 
state. For any two-qubit quantum state we can choose a 
set of 16 linearly independent operators {Mi} such that 
p can be decomposed as 



16 
i=l 



iMi. 



where the set {a} are the 16 parameters to be estimated. 
If the operators are observables, then the 16 expectation 



values nii — TT[Mip] determine cj by 

16 

rui = ^Tr[M,Mj]cj-. 

i=i 

Previous worki has shown that in cQED a homodyne 
measurement of the cavity is a faithful measurement of 
az ■ For a quantum bus with two qubits the measurement 
operator— is 



M = - f Q{t) dt 
Jq 



/3l(T^+/32< + /3l2Cr^ 



Here, Q is the measured quadrature amplitude, r is an 
averaging window, and the /3 are calibrated coefficients. 
For this experiment, r — 450ns and (/3i, /?2, /3i2) ~ 
(60,50,40) pY. 

Since the measurement contains both one- and two- 
qubit operators, a complete set of linearly independent 
operators Mi can be made by applying only single-qubit 
rotations prior to measurement. The set of 15 pre- 
rotations used in this experiment is all combinations of 
/, Rx, rV^ , Ry^^ on left and right qubits, except that 
Rx (8> Rx is not used. Only 15 measurements are needed 
to determine p because of the constraint of trace normal- 
ization, Tt p = 1 (equivalently we choose Mig — I, which 
always gives mie = 1). 

Experimental averages m; are obtained by recreat- 
ing the quantum state (executing the gate array), pre- 
rotating and measuring 450,000 times. While ideally p 
could be obtained from the experimental rUi by inver- 
sion of TilMiMj], this method pays no attention to the 
properties p must have: hermiticity and positive semi- 
definiteness (trace normalization is included by the choice 
of decomposition). These physical constraints are auto- 
matically included by a parametrization 

P ~ Tr[TtT] ' 

where T is a lower triangular matrixi^. For two qubits, 

( ti 0^ 

^5 + ih t2 

^11 + iti2 h + its h 

\tl5 + itiQ ti3 + iti4 tg + itiQ tiJ 



T 



2 



The ti are found by standard Maximum likelihood where the ai are weighting factors. We weight all mea- 
EstimationS with a likelihood function surements equally since amplifier noise dominates the er- 

ror in all the measurements. 

16 

£ = ^a,(m,-Tr[Af,p])^ 

i=l 
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Figure SI: Implementation of Deutsch— Jozsa algorithm, a. Gate sequence solving the Deutsch-Jozsa problem. The 
two-qubit gates Ui performing the transformation |Z,r) — > |Z,r © fi{l)) (© denotes addition modulo 2) for fo{x) = 0, fi{x) — 1, 
f2ix) = x, and faix) = 1-s are f7o = /®/, Ui = 7®i?J, U2 = {I®rV^RI)cUoo{I®RI^^), and Ui = {I®Ry'''^Rl)cUii{I® 
Ry^^^), respectively, b— e, Real part of the inferred density matrix of the algorithm output in the four cases (imaginary 
elements of Pmi are less than 0.05, 0.03, 0.05, 0.06, respectively). For the constant (balanced) functions fo and /i (/2 and 
fa), /5ml reveals high fidelity to |1,0) (|0, 0)), as expected. For the tomographs shown, the fidelity to the ideal output state is 
F = 0.94, 0.95, 0.92, and 0.85, respectively. Statistics for 8 runs of each of the four cases are given in Table I. 
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Figure S2: Experimental setup and wiring, a, Block diagram of room-temperature electronics. Arbitrary waveform 
generators, with 1 ns sampling rate and 10-bit resolution, produce voltages Vl and Vr, directly applied to the flux-bias lines, 
the I-Q modulation envelopes for the microwave tones driving single-qubit x- and y-rotations, and the pulse that modulates the 
cavity measurement. On the output side, an I-Q mixer and a two-channel averager (2 ns, 8-bit sampling) complete the readout 
chain performing homodyne detection of the cavity quadratures. The arbitrary waveform generators, microwave synthesizers 
and acquisition card are clocked with a Rubidium frequency standard (SRS FS725, not shown), b. Schematic of the microwave 
wiring inside the dilution refrigerator, showing heavily-attenuated input lines and an output chain with ~ 100 dB gain in the 
4-8 GHz range. 
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Figure S3: Microwave and flux pulses realizing the gates, a. An example sequence, executing the Grover search algorithm 
with oracle O — cUio and measuring Afis = — /^icr^ -I- /320"^ — /3i20"z (X" cr^. b. Illustration of the microwave and flux pulses 
realizing the operations directly above. All microwave pulses implementing the x- and y-rotations have Gaussian envelopes, 
with standard deviation a = 2 ns, truncated at ±2a. The rotation axis is set using I-Q (vector) modulation (see Fig. IS2[|. and 
rotation angle is controlled by pulse amplitude. A buffer of 5 ns is inserted between all microwave and flux pulses to avoid any 
overlap. 



